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We consider a supersymmetric model with extra U(1)b X U(1)l gauge symmetry that are broken 
spontaneously. Salient features of this model are that there are three different types of cold dark 
matter (CDM) candidates, and neutral scalar sector has a rich structure. Light CDM with crsi ~ 
1(T 3±1 pb can be easily accommodated by leptophobic dark matter (xb) with correct relic density, 
if the U(1)b gauge boson mass is around 2m XB . Also the PAMELA and Fermi/LAT data can be 
fit by leptophilic CDM with mass ~ 1 TeV. There could be interesting signatures of new fermions 
and new gauge bosons at the LHC. 
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Proton decay is possible only if both baryon and lep- 
ton numbers are broken. The current bound on proton 
life time indicates that both of these symmetries are ex- 
tremely good ones. This fact might indicate that these 
symmetries are exact local symmetry, rather than global 
symmetries which are always expected to be broken to 
some extent. 

Recently the gauged U(1)b X U(1)l model without su- 
persymmetry has been constructed and phenomenologi- 
cal consequences of this model have been discussed 
Anomaly can be cancelled by introducing new full fam- 
ily of the opposite (or the same) chiralities from the SM 
fermions with the baryon number — l(l)and lepton num- 
ber —3(3), respectively. Now, let us focus on the model 
with the opposite chirality. 

Basic features of the model with a CDM candidate 
are additional scalars, Sl and Sb that break the gauged 
U(1)l and U(1)b at some scale, and a new stable scalar 
X which makes new heavy quark decay into light quark 
and X. The model also has interesting aspects for neu- 
trino physics and leptogenesis []]]. 

The model considered in Ref. [H, being a nonsuper- 
symemtric model with fundamental scalars, does not ad- 
dress the fine tuning problem of Higgs mass 2 . If we con- 
sider the fine tuning problem seriously, we have to include 
another new physics at EW scale. Softly broken super- 
symmetry at EW scale is probably the simplest way to re- 
solve this issue. Furthermore, the lightest neutral particle 
can be stable in supersymmetric models with R-parity. 
The extension of the gauged U(1)b X U(1)l model to 
SUSY can realize not only R-parity spontaneously, but 
also additional CDM candidates. There is a possibility 
that we can provide several experimental results with a 
comprehensive explanation. 

In this letter, we consider SUSY extension of U(1)b x 
U(1)l model with new full family fermions with opposite 
chirality, and study its phcnomcnological consequences. 
(The model with the same chirality fermions is similar 
to the MSSM4, and we do not pursue this possibility in 
this letter.) The matter content of our model is similar 
to that of Ref. [l[, but we add lepton number carrying 
supcrfield Xl and Xl, which can make leptophilic CDM. 



Some of the salient features of our model arc that both 
the dark matter and neutral Higgs boson sector become 
very rich with multicomponent. There are three CDM 
candidates, leptophilic and leptophobic and the usual 
lightest neutralino. There are ten neutralinos and ten 
neutral scalars. After baryon and lepton numbers are 
broken spontaneously, there are mixing among new neu- 
tral fermions (scalars) and those in the minimal super- 
symmetric standard model (MSSM). 

In the SUSY U(1)b xU(1)l model, all gauge anomalies 
are cancelled by introducing one new family of quarks and 
leptons with opposite chiralities from the SM fermions. 
(See TablelH) Sb(Sb) and Sl(Sl), whose charges are 
ib(-Bb) and 2(— 2), are the fields which break U(1)b 
and U(1)l- In order to realize spontaneous U(1)b and 
U(1)l breaking, we set soft SUSY breaking terms of 
Sb(Sb) and Sl(Sl), squared mass terms and B-terms, to 
ones in appropriate parameter regions like MSSM higgs 
without assuming specific mediation mechanism. 

Sl and Sl couple with the right-handed neutrino, Ni, 
and generate heavy Majorana masses. Xb and Xb are 
the fields which make extra quarks, such as Q' , decay. 
If we assume that superpotential is a polynomial, extra 
quarks cannot decay without Xb and Xb- We forbid 
couplings between Sb(Sb) and SM particles according 
to the charge assignment in Table HI Compared with non 
SUSY model of Ref. [l[, we have additional Xl and Xl 
with U(1)l charges ±i%l which can be leptophilic CDM. 
In order to forbid Xl coupling to leptonic bilinear, we 
impose til 7^ 0, ±1, ±2, ± 4. The charge of Sb, 
tib, also satisfies ns ^ 0, ± 1/3, ± 2/3, ± 4/3 to 
forbid direct couplings between Sb{Sb) and Xb{Xb) in 
renormalizable superpotential. For nj, = ±1 and tib = 
±4/3, operators like SlX\ and SbX b are allowed and 
Xl and Xb could be stable, but we assume that ul 7^ ±1 
and ub ^ ±4/3 for convenience. 

Superpotential of the model is given by 



W = Wquark + Wle; 
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where the first two are Yukawa terms and the last one is 
the bilinear term: 
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TABLE I: SU{3) c xSU(2) l xU(1) y xU(1)bxU(1)l charges 
of fields in this model. i,j, k = 1,2,3 are the generation 
indexes for the SM fields, and the primed fields are newly 
added to cancel all the gauge anomalies. All fields are left 
chiral superfields. 



Field 


Charges 


Fields 


Charges 


Q l 


(3, 2, 1/6; 1/3,0) 


Q 

u' 


(3, 2, -1/6; -1,0) 


U* 


(3, 1,-2/3; -1/3,0) 


(3, 1,2/3; 1,0) 


D l 


(3,1, +1/3; -1/3,0) 


D 


(3, 1,-1/3; 1,0) 


V 


(1,2, -1/2; 0,1) 


L 


(1,2, 1/2; 0,-3) 


E l 


(1,1,1; 0,-1) 


E 


(1,1,-1; 0,3) 


N* 


(1,1,0; 0,-1) 


N 


(1,1,0; 0,3) 


H u 


(1,2, 1/2; 0,0) 


H d 


(1,2, -1/2; 0,0) 


X B 


(1,1,0; 2/3,0) 


As 


(1,1,0; -2/3,0) 


Sl 


(1,1,0; 0,2) 


sZ 


(1,1,0; 0,-2) 


S b 


(l,l,0;n s ,0) 


Sb~ 


(l,l,0;-n s ,0) 


X L 


(l,l,0;n L ,0) 


IZ 


(l,l,0;-ni,0) 



In this model, lepton and the baryon numbers are gauge 
charges, so that there are no usual R— parity violating 
couplings. Assuming that ub is 2k /3 (k ^ 0, ±i, ±1, ± 
2), where k is integer, U{1)b x U(1) l breaks to Z 2 b x Z 2 l- 
R-parity, (— 1) 3B + L + 2 J J is automatically conserved after 
U(1)b x U(1)l symmetry breaking. 

We assume that Xl and Xb do not get nonzero VEVs, 
in order to avoid large mixing between Q' and Qi which 
would generate large FCNC. In the lepton sector, the 
yukawa couplings, such as Xi, generate mass mixing be- 
cause of nonzero < Sl >■ For example, fi — > e, 7 process 
gives strong constraints on the yukawa couplings, 



< 



10 



-13 



(Sl)(Sl) GeV 



(5) 



There are two charged scalars as in the MSSM. The 
neutral scalar bosons are much richer than the MSSM: 

• CP-even : (Re (S B ) , Re (&)), (Re (S L ) , Re (3£)) 

• CP-odd: (Im(S B ),Im(SF)), (Im (S L ) , Im (ST)) 



After U(1)b X U(1)l is spontaneously broken, there are 
mixings between the MSSM neutral Higgs bosons and L 
and B charged neutral scalars (S B , S L and their conju- 
gates) at one-loop level. While MSSM Higgs bosons cou- 
ple to the fermion masses, the L or B charged scalars cou- 
ple to the baryon and lepton numbers of fermions. There- 
fore after making field redefinition to the physical mass 
cigcnstatcs, the neutral scalars can have decay patterns 
which are different from the MSSM Higgs bosons, espe- 
cially for light fermions for which m f / v sm. (3(cos P) ~ gs 
or t/£. The mixing and upper bound for the lightest 
U(1)b biggs will be very small, like O(10) GeV, in the 
parameter region we discuss below. We expect the light- 
est to evade constraints on colliders Q , but more detailed 
studies are in need 0. 

The neutralino sector of our model is also enlarged due 
to the new gauge symmetries. Let us call As and Xl the 
gauginos of U(1) B and U(1)l, respectively. There are 

two more sets of neutral fermions, (Xb, Sb, Sb) = ^b T i 
and (Xl, Sl, Sl) = , 3 / " T in addtion to the usual 3 neu- 
tralinos (B,W 3 ,H u ,Hd) = fgg (gauginos plus Higgsi- 
nos). All in all there are 10 neutral Majorana fermions. 
After the gauge symmetry of our model is spontaneously 
broken into SU(3)c X U(l) ein , all these neutral fermions 
mix with each other either at one-loop or two loop level. 
The mass eigenstates Xn=i 2 10 °f ^ en neutral fermions 
are linear combinations of ^ L and ^% H - 

Xl(Xl) and Xb(Xb) are also neutral, but they do not 
have mixing with Xb according to the zero VEVs of the 
bosonic superpartners. That is, they can be dark matters 
of Dirac fermion and respect accidental U(l) symmetry. 
Let us denote the dirac fermions as just Xb and Xl in 
the following sections. The bosonic sectors of Xl and Xb 
also have global U(l), so that they cannot have mixing 
with the other neutral scalars, or the mixing could be 
very small. 

Our model has multicomponent CDM's due to three 
conserved quantities: the usual R— parity and two acci- 
dental global symmetries, U(l)x B and U(l)x L - 



U(1) Xb : (X B ,X B ,Q ,U ,D 



(e la X B ,e-* a X B ,e- 



,e la U ,e ta D ) 



(6) 
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U(1) Xl : (X L ,X L ) -> (e i0 X L ,e 

I 



(7) 



The SM particles, extra quarks and leptons arc R— parity 
even, whereas their superpartners and all gauginos are 
R— parity odd. While the extra neutral bosons, such as 
Sb, are R- parity even, their superpartner is odd. Im- 
parity realizes a stable particle in our model, so that x?) 
Xb, and Xl are good candidates for CDM. 

Besides, both Xb and Xb carry U (l)x B charge Qx B = 
+1, and their conjugates carry Qx B = — 1- For lep- 
tonic fields, Xl,Xl (and their conjugates) carry Qx L = 
1(— 1), whereas all the other particles carry no U(l)x L 
charge. There are 2 complex scalars charged under only 
U(l)x B , but they are not degenerate generally because 
of B-term, BxXbXb- The heavier scalar field can decay 
to 2 quarks and the lightest scalar. Unless the light- 
est is heavier than both Xb and xl-* it cannot decay 
because of the global symmetry. The lightest U(l)x B 
charged scalar particle, Xbi, can be leptophobic CDM. 
The U(l)x L charged fields Xl\ (scalar) is also a good 
leptophilic CDM candidate. 

Eventually, there are 1 Majorana fermion, 2 Dirac 
fermion, and 2 complex scalar fields as good CDM candi- 
dates in this model. Among them, 3 particles can be sta- 
ble corresponding to 3 global symmetries, and the other 
2 fields decay to 2 CDMs. It depends on the mass spec- 
trum which fields are stable or not, and each sets must 
allow extra quarks to decay. The extra quark, such as 
Q' , is i?-parity even and charged under U(l)x B , so that 
at least lighter Xbi or JTg must be lighter than Q' . If 
Xbi is heavier, both of Xb and x? must be lighter than 
Q' because of i?-parity. 

Generally speaking, U(l)x B is anomalous and U(l)x B 



breaking terms, such as XbS 



-2/3n B 
B 



XbS b k (with 



|fc| 7^ 0,1/2,1,2), can be written down, including non- 
renormalizable terms^FIowever, such breaking terms al- 
low mixing between Xb and neutralino, because Xb gets 
nonzero VEV through tadpole terms like Xb(Sb)~~ ■ If 
our supcrpotcntial is a polynomial in supcrficlds, such 
dangerous operators would be forbidden as far as \k\ > 1 
is satisfied, and the stability of Xbi and dirac CDM, Xb, 
could be guaranteed. This argument can be applied to 
U(l)x L , and jil cannot be even integer and ±1 to respect 
U(l)x L . 

Since there are many new particles and parameters 
compared with the MSSM, it would be easy to imag- 
ine that there are rich phenomenology of this model 
both at colliders and in CDM. In order to illustrate the 
richness of our model, we choose to fit both (arguably) 
intriguing signature of light cold dark matter (1CDM) 
from direct detection experiments by CoGeNT Q and 
DAM A/LIBRA H, and positron excess observed by 
PAMELA @ and Fermi/LAT satellite 0). We consider 



the scenario that there are one light leptophobic CDM, 
Xb or Xbi, lightest neutralino, and one heavy leptophilic 
CDM, X L or X L1 . 

CoGeNT and DAMA/LIBRA results are challenged by 
CDMS and XENON10 @, so we need very care- 
ful analysis to realize them in our model [Io| . Even 
if they are excluded, there is a region close to the Co- 
GeNT and DAMA signal region, where 1CDM xb can 
have <7si ~ 10~ 3±1 pb with an acceptable thermal relic 
density, as long as Mz B /gB ~ ITeV and Mz B ~ 2m XB 
which is not excluded by collider data (see Fig. 1 and 
related discussions). 

Keeping this remark in mind, let us focus on the 
leptophobic 1CDM scenario around the CoGeNT and 
DAMA/LIBRA signal region. We assume that the 1CDM 
claimed by CoGeNT and/or DAMA/LIBRA is a lepto- 
phobic dark matter, Xb (Xbi), in our model. 

From the claimed cross section asi ~ 10 -40 cm 2 and 
the 1CDM mass ~ 7 GeV , we can obtain a constraint 
on a combination of mz B /gB ~ 1 TeV, if the contribu- 
tions of superparticles are negligible because of their TeV- 
scale masses. Using this constraint, we can constrain 
further the mass of mz B by calculating flicDAih 2 and 
assuming it is equal to or less than the WMAP measure- 
ment of ficDM- We find that the relic density turns out 
too large for heavy mz B ■ We can use the s— channel anni- 
hilation on the Zb resonance in order to enhance the an- 
nihilation cross section, thereby getting the right amount 
of relic density for the cold dark matter. Then there 
is more or less unique value for mz B ~ 2micDM ~ 14 
GeV, with very small U(1)b gauge coupling as ~ 10~ 5 . 
See Fig. 1 for detail. The X B i CDM does not have S- 
wave contributions, so that the allowed region is nar- 
rower. Squark exchanging could also contribute to the 
relic density through the t— channel annihilation, but it 
would be difficult to enhance the annihilation cross sec- 
tion. 

Such a light U(1)b gauge boson with a weak gauge 
coupling is definitely out of reach at the current searches 
from colliders and low energy processes 12j, [13( . It would 
be interesting and important to search for such light lep- 
tophobic gauge boson at the upcoming experiments. 

Positron excess observed by PAMELA collaboration Q 
might be due to dark matter pair annihilation into 
charged lepton pairs. There are explicit models based on 
leptophilic U(l) symmetries 14lll5|. Since our model has 
lepton charged carrying neutral scalar (Xl\) or fermion 
(Xl) that couples only to leptons due to the gauged 
U(l) l and supersymmetry, we could try to fit PAMELA 
data using this leptophilic dark matter. 

In order to realize PAMELA, we have to consider large 
annihilation cross section, compared with (av) <~ 1 pb 
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FIG. 1: Allowed region in the (log ob, Mz B ) plane for 
miCDM — 7 GeV, that is consistent with WMAP data, is 
the inside of each 2 line, where the leptophobic 1CDM is 
bosonic Xbi (dashed lines) or fermionic A'g (solid lines). 
Inside of thick black line is excluded by hadronic decay of 
T(liS'). The bands within two blue lines and two gray lines 
are DAM A/LIBRA and CoGeNT regions [TJ. 

which is favored by WMAP data, VLh 2 ~ 0.11. In our 
model, its pair annhilation channels will be 

X L X{ -4 Z* L -> ll.ox Z L Z L — » 41' s. 

The latter is dominant if Zl is lighter than Xl- In 
that case, Sommerfeld enhancement could work at low 
energy like (av{T )) ~ S(av(T )) . In the Fig. gj we 
can see the enhancement and the relic density with 
tbcdm = 3TcV. However, and Mz L are strongly con- 
strained by the experiments: e, scattering constrains 
MzJgL as M ZL /g L > 1.4 TeV 0, and the LEP-K 
bound, Mz L /gL ^ 6 TeV, which is discussed in U(1)b-l, 
could be applied to our model [l7j |. 

We notice that it is difficult to find the favored re- 
gion with large S according to Fig. [2] The heavier Xl 
with large S could avoid the constraints, but such heavy 
Xl pair annihilations may produce high energy neutri- 
nos (upto the mass of parent Xl) of 3 flavors at the same 
rateSjWhich may cause conflict with the cosmic observa- 
tion [lBj]. 

Another CDM also contributes to the relic density, but 
if it has a heavy mass, like lOOGcV, the contribution 
could be enough small, as several works have discussed 
so far. For example, it has been discussed that yukawa 
coupling, XbXBib'b, can realize not only small relic den- 
sity of Xbi but also the favored direct scattering cross 
section in Ref. 0. If Xi is MSSM-like CDM, Z boson 
and Higgs exchanging are helpful . 

Kinematic mixing terms among U(1)b, U(1)l and 
U(1)y are allowed by those gauge symmetries, and gencr- 
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FIG. 2: This figure shows the relic density and Sommerfelt 
enhancement of Xl with menu = 3 TeV and Qx L = 1. 
The above of the blue line corresponds to (yih?)x L < 0.11. 
In the white, yellow and gray regions, the enhancements are 
S > 100, 100 > 5" > 10, and 10 > S. The left sides of the 
dashed red and the red lines are exculded by e, f M scattering 
and LEP-H bound. 

ated radiatively. Their upper bounds are very tight [l3| . 
but the radiative correction is very tiny if gauge cou- 
plings are small enough as in our scenario. In any case, 
the initial condition must be controlled to suppress the 
kinematic mixing, as discussed in Ref. [l3| . 

The conserved quantities are U(l)x B > U(l)x L charges 
and R— parity. Note that Xb and Xl are R— parity even, 
whereas their fermionic superpartners are R— parity odd. 
Let us discuss some decays of new particles as examples 
(Xi is the lightest R-): 

X B -> qq* -> q(qXi) (8) 
X B -> qqXiXs (9) 

Sb and its conjugate mix with Xb and decay intoj7g* — > 
q(qXi)- Similarly, Al, Sl andXt decay into UxiXl- Un- 
like the MSSM, the R— parity even scalar Xb can decay 
into qqxiXs through two processes: 

X B -> ^x B X B -> X B {qqXi) 
-> Q'q* -> (qXB)(q jXl ) 

if kinematically allowed. However Xb can be stable if 
this decay is kinematically forbidden. This is guaran- 
teed by conservations of U(l)x B , U(l)x L charges and 
R— parity. If Xb is heavier than Xb, the reverse pro- 
cess will occur if kinematically allowed. The new heavy 

(s)quarks can be produced at the LHC: pp — s- 1 t' , t't' . If 
Xb is lighter than Xb, the main decay modes of t and 
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t' are t qX and t' — > Xgq — > Xgqxi, respectively. If 
Xb is lighter than Xb, then the main decay modes of t 
and t' arc t — > Xsq — > XsqXi, an d ^ ~~ ^ ^bI- If some 
of the decays are kinematically forbidden, t or t' will be 
stable and hadronized, and thus strongly constrained by 
heavy charged particle search. Finally, new scalars Sb, 
Sl and their conjugates will mix with neutral Higgs and 
will decay into lighter fermions or weak gauge bosons or 
Z B 's. 

There are strong constraints on t and b masses. They 
can decay to quarks and CDM, but such exotic quarks 
search could be replaced by SUSY particle ones, such 
as squark decay process, t — > t. In Ref. [2fJ, the 
bounds for t and b masses, which correspond to the 
case with scalar Xb CDM, are dicussed in that way: the 
lower masses are more than 300GeV. If this argument is 
applied to our model, Landau poles of Yukawa couplings 
appear at TeV scale. 

We must also care about EWPT bound for higgs mass, 
extra quark masses, and all sparparticlc masses. The 
discussions in Ref. [2l[ can be applied if superparticle 
and extra higgs contributions are small. 

In this letter, we constructed supersymmetric exten- 
sion of U(1)b x U(1)l model, extending the works by 
Wise et al. [1|. The model is anomaly free, with very rich 
phenomenology in terms of colliders, Higgs bosons, flavor 
physics and CP violation, and most notably CDM. The 
most interesting aspect of this model is the multi com- 
ponent CDM's. The model houses both leptophilic and 
leptophobic CDM's in addition to the usual lightest neu- 
tralino, and the stability of these CDM's are consequence 
of gauge symmetry and particle contents, not of ad hoc 
Z2 symmetry. Furthermore, the leptophobic CDM can 
be light enough to accommodate the CoGeNT/DAMA 
signal region, with correct amount of relic density. The 
leptophilic CDM can be heavy enough to explain the 
positron excess observed by PAMELA and Fermi/LAT. 
However, it seems to be difficult to explain PAMELA 
without assumption of large boost factor. 

We could describe the simpler scenario that the lightest 
neutralino approximately becomes leptophilic because of 
small mixing, and realizes PAMELA data. Furthermore, 
our model has an interesting signature that the missing 
Et signature could be due to two different DM (with 
different masses and spins). We touched relevant features 
of our model with a wide brush, relegating the details of 
phenomenology to the future works [3[ . 
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